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ABSTRACT 

Aims. The properties of dust attenuation at rest-frame UV wavelengths are inferred from very high-quality FORS 2 spectra of 78 
galaxies from the GMASS survey at redshift 1 < z < 2.5. These objects complement a previously investigated sample of 108 
UV-luminous, intermediate-mass (i.e., with stellar masses ~ 10'" - 10" Mq) galaxies at similar redshifts, selected from the FDF 
spectroscopic survey, the K20 survey, and the GDDS. Detection of the broad absorption feature centred on about 2175 A ("UV 
bump") implies that the average UV extinction curve of a galaxy more closely resembles that of the Milky Way (MW) or the Large 
Magellanic Cloud (LMC), and differs from that of the Small Magellanic Cloud (SMC). 

Methods. The shape of the UV extinction curve is constrained by a parametric description of the rest-frame UV continuum with the 
support of a suite of models combining radiative transfer and stellar population synthesis. The UV bump is further characterised by 
fitting Lorentzian-like profiles. 

Results. Spectra exhibit a significant 2175 A feature in at least 30% of the cases, especially those suffering from substantial reddening. 
If attenuation is dominated by dust ejected from the galaxy main body via galactic winds or more localised superwinds, UV extinction 
curves in-between those of the SMC and LMC characterise UV-luminous, intermediate-mass galaxies at 1 < z < 2.5. The fraction of 
galaxies with extinction curves differing from the SMC one increases, if more dust resides in the galactic plane or dust attenuation 
depends on stellar age. On average, the width of the manifested UV bumps is about 60% of the values typically measured along 
sightlines in the LMC and MW. This suggests the presence of dust similar to that found in the "supergiant" shell of ionised filaments 
LMC 2, close to 30 Dor. The presence of the carriers of the UV bump (probably organic carbon and amorphous silicates) at 1 < z < 2.5 
argues for outflows from asymptotic giant branch (AGB) stars being copious then. Consistent with their higher star-formation rates, the 
GMASS galaxies with a manifested UV bump are more luminous at rest-frame 8 yum, where the emission is dominated by polycyclic 
aromatic hydrocarbons (also products of AGB stars). In addition, they exhibit larger equivalent widths for prominent UV (metal) 
absorption features, mostly of interstellar origin, which indicates overall more evolved stellar populations. 

Conclusions. We conclude that diversification of the small-size dust component (responsible for the differential extinction at UV 
wavelengths and the emission at mid-IR wavelengths) has already started in the most evolved star-forming systems at 1 < z < 2.5. 

Key words, galaxies: high-redshift - galaxies: starburst - galaxies: ISM - ISM: dust, extinction - ultraviolet: galaxies 



1. Introduction 

A broad absorption excess centred on 2175 A (the so-called "UV 
bump") characterises the extinction curves usually inferred for 
sightlines towards the Milky Way (MW) and Large Magellanic 
Cloud (LMC) (for reviews see Fitzpatrick 2004 and Clayton 
I2004I I. and also towards M31 (Bianchi et al. |19961 l. In contrast, 
this feature is absent in the typical extinction curves of the Small 
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Magellanic Cloud (SMC), which exhibit a steep rise at far-UV 
wavelengths. SMC-like dust appears to characterise starburst 
galaxies in the local Universe (Calzetti et al. 119941 ) and Lyman 
break galaxies (LBGs) at redshifts 2 < z < 4 (Vijh et aL |^03ll 
For this reason, the so-called "Calzetti law" (Calzetti et al. 11 9941 
2000) is frequently assumed to describe attenuation by dust in 
starburst and star-forming galaxies in the distant Universe. 

However, there is strong support for the presence of the 
2175 A feature in spectra of at least 30% of the vigorously star- 



2 



S. Noll et al.: GMASS ultradeep spectroscopy of galaxies at z ~ 2. IV. The variety of dust populations 



forming galaxies at 2 < z < 2.5 (Noll & Pierini 2005, here- 
after referred to as NP05) and 1 < z < 1.5 (Noll et al. 2007, 
hereafter referred to as NPPS). Furthermore, an LMC-like ex- 
tinction curve with a UV bump was inferred clearly from both 
the spectrum and the broad-band spectral energy distribution 
(SED) of the afterglow of the gamma ray b urst G RB 070802 at 
z = 2.45 (Fynbo et al. 2007; Kruhler et al.'2008'i Eliasdottir et 
al. [2008 ). As for other objects at intermediate/high redshifts, the 
presence of the carriers of the 2175 A feature is controversial 
in quasars/active galactic nuclei (AGN) (cf. Pitman et al. 120001 
Maiolino et al. lIOUTl Hopkins et al.[20 Ml Ver net et al.[200 T]) and 
intervening Mgll absorbers (Malhotra [19921 Wang et al. 120041 
Vanden Berk et al. 2008) , whereas it is likely in Ca II absorbers 
at intermediate redshifts (Wild & Hewett l2055T l. 

It would therefore be important and interesting to under- 
stand when the carriers of the UV bump first appeared in the 
Universe and how they can survive in vigorously star-forming 
galaxies at high redshifts. Multiple carriers (i.e., organic carbon 
and amorphous silicates) are likely to explain the almost invari- 
ant central wavelength and variable bandwidth of the 2175 A 
feature (Bradley et al. 2005 and references therein). However, 
a random, hydrogen-free assembly of microscopic sp^ carbon 
chips, forming a macroscopically homogeneous and isotropic 
solid, has recently been proposed as a model carrier of the UV 
bump by Papoular & Papoular ( I2009I I. According to this promis- 
ing model, different grain formation histories and/or different in- 
terstellar radiation fields and other heating effects are responsible 
for a change in bandwidth without notably impacting the central 
wavelength. Whatever the carrier of the UV bump, the propor- 
tions of the interstellar dust ingredients that are responsible for 
both the UV bump and the far-UV rise in the extinction curve, 
relative to each other, are approximately similar in both the MW 
and LMC, but appear to differ significantly in the SMC (Whittet 
12003) . Differences also exist within the MW for different sight- 
lines and towards different environments, from dense molecular 
clouds to the diffuse interstellar medium (ISM; see Fitzpatrick 
[2004). They must reflect the extreme sensitivity of the small-size 
dust components to the local chemical enrichment and energy 
budget (in terms of radiation field and shocks), as well as their 
selective removal from the size distribution due to a number of 
physical processes (e.g., Whittet i2003t Gordon et al. 12003b . 

We continue an analysis started by INP05I and extended by 
INPPSI We adopt a new homogeneous sample of UV-luminous, 
intermediate-mass galaxies at 1 < z < 2.5 drawn from the 
spectroscopic catalogue of the galaxy mass assembly ultra-deep 
spectroscopic survey (GMASS, Kurk et al. |2p08a; Cimatti et 
al. 2008; Kurk et al., in p rep.) i n the GOODS-South Field 
(GOODS-S; Giavalisco et al. |2UU4b. T his sample has about 70% 
of the size of the entire sample in NPPS but contains mostly ob- 
jects in the poorly investigated range 1 .5 < z < 2.0, known as the 
"redshift desert". The excellent dataset available for the GMASS 
sample enables us to characterise more accurately the appar- 
ent shape of the 2175 A feature in high-redshift star-forming 
galaxies, and to investigate the relation between dust properties 
inferred from attenuation at UV wavelengths and emission at 
mid-IR wavelengths. Furthermore, it adds relevant information 
about the link between small-size dust populations (i.e., large 
molecules and small grains) and chemical abundances in stars 
and gas of the host galaxies. 

Throughout this paper, we adopt a ACDM cosmological 
model where Qa = 0.7, Om - 0.3, and Hq - 70kms ' Mpc 
By default, photometry is given in the Vega magnitude system; 
otherwise AB magnitudes are indicated. 



2. The sample 

2.1. The GMASS sample 

The GMASS spectroscopic catalogue (Kurk et al. |2008al Kurk 
et al., in prep.) comprises 167 objects with secure spectroscopic 
redshifts. It is based on a mass-sensitive selection criterion using 
photometry in the Spitzer 4.5 ;um IRAC filter (mp^^iA.S fim) < 
23), and the additional photometric -redshift criterion Zphot > 1 -4. 
For our analysis, we identify and select vigorously star-forming 
objects, for which we can robustly describe the continuum emis- 
sion at rest-frame wavelengths from 1700 to 2600 A in the para- 
metric way illustrated bv lNP05l and lNPPSI (see Sect.ITT]). In the 
(spectroscopic) redshift range of our interest (1 < z < 2.5), 
GMASS yields a total of 78 class 2 objects, i.e., actively star- 
forming galaxies with UV spectra dominated by young stellar 
populations and no AGN features (see Mignoli et al. 2005 ). For 
68 objects spectroscopically confirmed to lie at 1.5 < z < 2.5, 
selection biases are understood: the average galaxy is about 0.2- 
0.25 mag fainter at 4.5 yum ((mAB(4.5 /vm)) = 22.0, cr - 0.7) than 
the average objects in the full spectroscopic and photometric 
samples for ;«ab(4.5 ^m) < 23 and 1.5 < z < 2.5. Conversely, it 
is about 0.3 and 0.8 mag brighter in the visual ((Veofe.AB) = 24.3, 
cr - 0.5) than the average galaxies in these two catalogues, re- 
spectively. Therefore, the following analysis addresses almost 
exclusively the sample of 68 GMASS UV-luminous galaxies at 
1.5 <z<2.5. 

The spectroscopic observations for the GMASS sample were 
carried out with FORS 2 at the VLT using the blue 300V and the 
red 3001 grisms, at resolutions of 440 and 660, respectively. Only 
a minority of the objects was observed in both masks. In view of 
the particularly long exposure times per mask of 1 1 to 32 h (and 
consequently high signal-to-noise ratio - S/N), the main reason 
for rejecting a spectrum is the limited wavelength coverage. 

For the GOODS-S, a wealth of multiwavelength data is 
available as part of the GOODS legacy (Dickinson et al. 120031 
Giavalisco et al. l2004l l. In particular, the deep Spitzer MIPS ob- 
servations at 24/im (Chary et al., in prep.) and IQjjm (Frayer 
et al., in prep.) are important to probe re-emission by dust at 
rest-frame mid-IR wavelengths. This analysis is also based on 
the available broad-band photometry from the UV to the near- 
IR (Kurk et al., in prep.), since it utilises galaxy parameters 
derived from model fits to the photometric data. For instance, 
both the star-formation rate (SFR) and the total stellar mass 
(M*) of a galaxy used here were derived from fitting broad- 
band SEDs (extending from U band to the IRAC 8 /im band) 
using the Hyperz code of Bolzonella et al. (I2000I I. the redshift 
being fixed to the measured spectroscopic value (cf. Pozzetti et 
al., in prep.). In particular, fitting made use of the stellar popula- 
tion synthesis models of Maraston ( |2005n including a Salpeter 
(I1955I I initial mass function (IMF) and fixed solar metallicity 
(Z0), whereas age was allowed to vary up to the lookback time 
of the object. Star-formation histories (SFHs) were described by 
exponentials with e-folding times between lOOMyr and infinity 
(constant SFR). Dust attenuation was assumed to be described 
as in Calzetti et al. (I2000I I. 

Finally, we note that morphologies of the GMASS galaxies 
were determined by a Galfit (Peng et al. ,2002) analysis of the 
HST-ACS imaging at 8500 A from the GOODS-S (Cassata et 
al., in prep.). For 1.5 < z < 2.5, the z-band maps rest-frame 
mid-UV to near-UV wavelengths, so that the morphological pa- 

' We use these particular models because they are the only published 
models including the thermally pulsating asymptotic giant branch phase 
for intermediate-mass (2-5 Mq) stars. 
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Fig. 1. R-band magnitude versus redshift for the total sample 
of 186 actively star-forming galaxies at 1 < z < 2.5 selected 
from the GMASS (filled squares), FDF spectroscopic survey 
(lozenges and crosses), the K20 survey (circles and plus signs), 
and the GDDS (asterisks). Crosses, plus signs, and asterisks 
mark galaxies without a determination of the UV continuum 
slope. To illustrate selection biases, the curves show lines of con- 
stant luminosity for typical 1 < z < 1 .5 FDF (solid line) and 
K20 (dashed hne) spectra. Both curves intersect at /? = 23.5 and 
z= 1.5. 

rameters reflect the distribution of young stellar populations. In 
this study, we adopt measurements of the effective radius r^, the 
Sersic (1968) index Wser, and the axial ratio b/a. Moreover, we 
investigate the CAS parameters concentration C and asymmetry 
A measured according to a method described by Cassata et al. 
(I2005l l. 

2.2. The FDF, K20, and GDDS samples 

The GMASS sample is complemented by the spectroscopic 
sample investigated by ,NPPS. consisting of 108 actively star- 
forming galaxies at 1 < z < 2.5. Of these 108 objects, 66 galax- 
ies were selected from the (mostly) /-limited spectroscopic cata- 
logue of the FORS deep field (FDF) spectroscopic survey (Noll 
et al. 120041) and originally inspected for z > 2 by NP05 , In ad- 
dition, 34 galaxies were taken from the A's < 20-selected K20 
survey in the GOODS-S and a field around quasar 0055-2659 
(Cimatti et al. 120021 MignoU et al. 120051 1. and 8 galaxies were 
selected from the /- and /Tj -limited Gemini deep deep survey 
(GDDS, Abraham et al. 2004). For 88 objects from the FDF and 
K20 samples, the UV continuum slope could be determined in 
addition to the presence/absence of the UV bump (see lNPPST ). 

2.3. Basic sample properties 

Figure [T] shows the range of R-band magnitudes and redshifts 
covered by the different samples considered in this study, includ- 
ing the 10 GMASS UV-luminous galaxies at 1 < z < 1 .5. All but 
two lNPPS objects are brighter than R - 24 for z < 2, and a few 
reach R - 24.7 at higher redshifts. With GMASS, fainter galax- 
ies (24 < R < 25) represent half of the sample ((R) = 23.9) 
and are probed for the entire redshift range studied. Additional 
differences in the redshift distributions exist, as previously an- 
ticipated: 76% of the K20 galaxies are at 1 < z < 1.5, whereas 
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Fig. 2. R - Ks colour versus /Tj magnitude for sample galaxies 
with 1.5 < z < 2.5. Symbols are the same as in Fig.[T] 



the FDF sample mainly consists of objects at 1 < z < 1.5 (38%) 
or 2 < z < 2.5 (52%); in conti-ast, 67% of the GMASS objects 
are located at 1 .6 < z < 2. 1 . As indicated by curves of fixed lu- 
minosity in Fig.[T] the brightness requirements for spectroscopy 
(and possible evolution in the galaxy population) cause an un- 
avoidable increase in the average rest-frame UV luminosity of 
the sample with redshift. Therefore, redshift-dependent results 
have to be interpreted with care. 

As additional information. Fig. |2] illustrates the R - K^^-Ks 
colour-magnitude diagram of the sample galaxies at 1 .5 < z < 
2.5, where most of the selected GMASS galaxies lie. There is 
good overlap between GMASS and NPPS objects at these red- 
shifts, the GMASS galaxies being only 0.2 mag brighter and 
0.3 mag redder than the coeval FDF galaxies (Ks - 20.8 versus 
21.0, and R - - 3.1 versus 2.8, respectively). Hence, in spite 
of the different selection criteria (i.e., based on the rest-frame 
near-IR and UV domains for the GMASS and FDF samples, re- 
spectively), the photometric properties of the selected GMASS 
galaxies are similar to those of the selected NPPS galaxies. As 
Fig. |2] demonstrates, no extremely red object (e.g., Pozzetti & 
Mannucci l20001 ) is present in the full sample (cf. Fig.|2]). Aspects 
of dust attenuation in these objects were discussed by Pierini 
et al. (I2004a) and Pierini et al. (120051 ). From the results of 
these analyses, we conclude that our sample galaxies are proba- 
bly neither heavily reddened, young starbursts nor normal star- 
forming, bulge-dominated systems. In principle, the lack of faint 
red galaxies in our sample due to brightness limits could affect 
the results of our analysis. However, investigations of the rela- 
tion between mass and extinction (Greggio et al. 120081 ) suggest 
that relatively massive star-forming galaxies with such proper- 
ties could be rare by nature. Hence, we can reasonably assume 
that not detecting such galaxies (if any) does not bias our con- 
clusions. 



3. Data analysis 

In this section, we briefly describe our method for constraining 
the shape of the UV extinction curve from optical spectroscopy 
of high-redshift galaxies. In a first step, constraints are obtained 
from a suitable parametric description of the rest-frame UV SED 
of a galaxy (Sect. 13. It . as extensively discussed in INP05I and 
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Fig. 3. The UV-bump indicator 734 versus the UV redden- 
ing measure fi\, for the GMASS (filled squares), FDF (open 
lozenges), and K20 galaxies (open circles) at 1.5 < z < 2.5, 
and a comparison sample of 24 local starburst galaxies (crosses; 
seeNP05). The diagram also shows different dust attenuation 
models (see legend) for a Maraston (I2005i l model with Salpeter 
IMF, continuous SF, an age of 100 Myr, and solar metallicity. 
The symbols are plotted in intervals of AEg-v - 0. 1 . A Maraston 
model with a SF e-folding time of 3 Gyr and an age of 3 Gyr 
would produce intrinsically redder continua and, thus, would 
slightly shift the plotted model curves by 0.20 to the right and 
by 0.08 to the top. 



INPPSI In a second, new step, the detailed shape of the observed 
UV bump is derived by fitting a model (Sect. l3.2l i. 

3.1. Parametric description of the UV continuum 



INP05I introduced a parametrisation of the rest-frame UV SED 
of a star-forming, dusty galaxy, which highlights the presence of 
the UV bump, and offers constraints on the far-UV slope of the 
extinction curve. It is based on power-law fits to different sub- 
regions of the UV continuum, in a similar way to the parametri- 
sation used by Calzetti et al. (|1994| l for measuring the UV con- 
tinuum slope p. Hereafter, we utilise only the parameters 734 and 
j6b (introduced by NPPS). The first parameter represents the dif- 
ference between the power-law slopes 73 and 74 measured in 
the wavelength ranges 1900-2175 A and 2 175-2500 A, respec- 
tively, excluding strong narrow features. By construction, 734 is 
particularly suitable for identifying the presence of the broad ab- 
sorption feature centred on about 2175 A: the more negative 734, 
the higher is the probability that a UV bump of some strength 
is present (see Fig. [3]). Moreover, 734 depends only very little on 
stellar population properties (NP05). Negative values of 734 are 
almost impossible for combinations of stellar populations atten- 
uated by bump-less reddening laws only. INP05] chose 734 = -2 
as a fiducial threshold for a robust detection of the UV bump (see 
INPPSI for an updated discussion of this threshold). The second 
parameter ySb, measured at 1750-2600 A, indicates the contin- 
uum slope across the UV bump, and is therefore a reasonable 



proxy for the dust reddening of the UV continuum, as well as 
P (cf. lNPPSK As in 'NPPS ', /3b is preferred to p (measured down 
to Lya) because the latter cannot be derived from optical spectra 
for galaxies at z < 2. By using fib for all galaxies across the entire 
redshift range 1 .0-2.5, we are able to compare their properties in 
a consistent way. 

We know that 734 alone does not infer whether the extinction 
curve is similar to either the MW or LMC curve, because it in- 
dicates the apparent amplitude/area of the UV bump, which is 
a function of the dust opacity (see, e.g., Witt & Gordon 2000i 
Since the dust reddening of UV SEDs depends on the opacity 
as well, it is the distribution of a galaxy in the fib-Ju plane that 
identifies the extinction curve (see Fig.[3]l. 

3.2. Lorentzian-like fitting of the UV bump 

The parameter 734 can easily be measured in spectra with a rel- 
atively narrow spectral coverage, and even in noisy ones, but 
it loosely constrains the observed shape of the UV bump, and, 
thus, the kind of extinction curve (i.e., MW-, LMC-, or SMC- 
like). The high S/N and wavelength coverage of the GMASS 
spectra offer the opportunity for a more accurate description of 
the UV bump. In particular, for z ~ 2 the UV bump sits almost 
in the middle of the observed wavelength range, which allows 
the bump to be fitted in its shallowest wings. For optimum re- 
sults, we investigate stacked spectra (introducing subsampling) 
and not individual ones (see Sect. |4.2| i. 

According to Fitzpatrick & Massa (1990' FM), the shape of 
the UV bump present in extinction curves is most accurately de- 
scribed by a Lorentzian-like "Drude" profile 



D{x, X(i, 7) = 



(x^ - xo ) + x^y 



„2-,,2 ' 



(1) 



where x, xq, and 7 are expressed in units of inverse wavelength, 
X() being the central wavelength and 7 the inverse of the FWHM 
of the UV bump. In the FM parametrisation, extinction curves 
k{A-V) are normalised by the reddening Eb-v, and the UV bump 
is reproduced by C3 D{x, xq, 7), where C3 is an amplitude (in units 
of inverse wavelength in quadrature). In later work, Fitzpatrick 
& Massa (I2007I I defined two additional parameters, based on 
combinations of the basic parameters xq, 7, and C3. In particular, 
£bump = C3/7^ is the maximum height above the linear baseline 
extinction, and Abump = ^cj/il-y) indicates the area subtended 
by the UV bump. 

In analogy with the FM formalism, we characterise the ob- 
served UV-bump in galaxy spectra as an excess with respect to 
a baseline attenuation. At variance with extinction curves mea- 
sured for sightlines to individual stars, the observed UV bump in 
a galaxy spectrum is affected by radiative transfer effects since 
spectra probe large regions of galaxies (see, e.g., Ferrara et al. 
1999; Pierini et al. |2004bl in case of photometry). We set the 
baseline to be the attenuation curve of nearby starburst galaxies 
(Calzetti et al. I2000I I. where the wavelength dependence is fixed 
and the amount of attenuation at a given wavelength scales with 
the observed reddening Eb-v at optical wavelengths. Since the 
curvature of the selected reference curve for the continuum at- 
tenuation is not supposed to differ strongly from the unknown 
real shape within the range of the UV bump, the choice of 
the baseline mostly affects the measured amplitudes of the UV 
bump, while centre and width are expected to be almost invari- 
able. Our approach is unorthodox but allows us to link our results 
to those obtained from studies of galaxy evolution, which adopt 
the Calzetti law as a description of dust attenuation in galaxies 
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at cosmological distances. Therefore, we refer to the parameters 
derived for our sample galaxies as "FM-like". 

In our procedure, regions of the spectra that exclude the 
2175 A featur^ and other known strong spectral lines are fitted 
by different Maraston (2005 ) models with exponentially decay- 
ing star-formation rates and reddened by a Calzetti et al. ( I2000I I 
attenuation law. By doing so, we are able to interpolate across 
the UV bump. Taking into account uncertainties and approxi- 
mations, fits are considered to be "good" if their significance 
reaches at least 95% of that of the best-fit model. For these good 
fits, the residual 2175 A feature (if any) is detectable in a spec- 
trum created by obtaining the ratio of the initial spectrum to the 
model. Retrieving the properties of the observed UV bump there- 
fore depends only on the best-fit Eb-v value of the Calzetti law. 
The uncertainty in Eb-v mainly affects the observed UV-bump 
amplitude C3, whereas xq and 7 are less dependent on the best-fit 
model. These FM-like parameters are derived from Lorentzian 
fits to each continuum-subtracted UV bump. Lorentzian fits are 
considered as "good" if their significance is at least 95% of that 
of the best-fit model. They are used to compute the mean and 
standard deviation of each FM-like parameter. To estimate the 
influence of noise on these fits, a set of synthetic spectra modi- 
fied on the basis of the error function are computed and analysed 
in the same way as the real data. The final variance in the FM- 
like parameters measured for a stacked spectrum is given by the 
standard deviation of the good fits and the statistical noise, added 
in quadrature. 

4. Results 

We first present and discuss results obtained from the parametri- 
sation and modelling of the rest-frame UV continua of our high- 
redshift galaxies (Sect. 14.11 ). whereas detailed shapes of the ob- 
served UV bumps are illustrated in Sect. 14.21 The link between 
carriers of the UV bump and properties of the parent galaxies 
is addressed in Sect. 14.31 Finally, the existence of a relation 
between the presence of these carriers and the amount of re- 
emission by dust at rest-frame 8 fim is discussed in Sect. 14.41 

4.1. Extinction curves witli a UV bump are not unusual 
among liigli-redsliift gaiaxies witli gas and dust 

Figure |4] illustrates the redshift distribution of the proxy for 
the presence of a UV bump in spectra, 734, and the contin- 
uum slope in the UV-bump region, /3\,, for the GMASS, FDF, 
K20, and GDDS samples. A total of 186 actively star-forming, 
intermediate-mass galaxies at 1 < z < 2.5 is shown, including 
the 10 GMASS objects at 1 < z < 1.5. As from INPPSI the 
fraction of UV-luminous galaxies at these redshifts that clearly 
exhibit a 2175 A feature in their spectra (i.e., 734 < -2) is about 
30%. In particular, this fraction is 25% for the combined FDF 
and K20 samples at 1 < z < 1.5, and 36% for the FDF sam- 
ple at 2.2 < z < 2.5. A similar fraction of 24% significant 
UY bump galaxies (SUBGs) characterises the GMASS sample 
at 1.5 < z < 2.5, 13 out of 17 detections being for objects at 
1 .5 < z < 2.2. This fraction drops to 10% for the GMASS galax- 
ies at 1 < z < 1.5. 

We verified that, for the GMASS sample (see Sect.|2]l, sub- 
sampling in terms of rest-frame UV luminosity (i.e., splitting 
the sample at R = 24.0) or large scale structure (see the over- 
density of galaxies at z ~ 1.61; Kurk et al. I2008bl ) does not 
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Fig. 4. The UV-bump indicator 734 (upper panel) and the UV 
reddening measure /3b (lower panel) versus redshift for 78 
GMASS galaxies at 1 < z < 2.5 (filled squares) in compari- 
son to 108 galaxies from INPPSl (20 without /3h measurement). 
Symbols are the same as in Fig.[T| 




Fig. 5. Average UV reddening measure /3b for unit bins of the 
UV-bump indicator 734. The following subsamples are shown: 
GMASS for 1.5 < z < 2.5 (squares and soHd fines), FDF for 
1.5 < z < 2.5 (lozenges and solid lines), FDF for 1 < z < 1.5 
(lozenges and dashed lines), and K20 for 1 < z < 1.5 (circles and 
dashed lines). Only bins with at least two data points are plot- 
ted. The jSb uncertainties shown are mean errors derived from 
the sample variance and the uncertainties of individual measure- 
ments. 



^ A region across the UV bump of width 900-1000 A is masked to 
avoid even faint wings. 
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impact this fraction. Hence, we conclude that the non-detection 
of SUBGs at 1.5 < z < 2.5 in the FDF and K20 samples 
(see NPPS ) is probably due to low number statistics, plus pos- 
sible selection and detection biases. In support of this conclu- 
sion, we note that the distributions of the GMASS, FDF, and 
K20 galaxies at 1.5 < z < 2.5 in the fih-y^A diagnostic dia- 
gram are similar (Fig. [3]). In particular, the average value of 734 
for the GMASS galaxies at 1.5 < z < 2.5 is consistent with 
that of coeval FDF galaxies (-1.01 + 0.17 versus -0.70 + 0.30) 
within the errors. The average value of the reddening in the 
UV, y8b, is slightly higher (-1.48 + 0.05 versus -1.71 + 0.07). 
At 1 < z < 1.5, the combined FDF and K20 samples exhibit 
stronger 2175 A features ((734) = -1.23 + 0.23) and redder UV 
continua ((/8b) - -1.23 ± 0.08) due to a population of galaxies 
with extreme UV-continuum properties, which is not present in 
our sample at higher redshifts (see Fig.|4]i. 

The GMASS sample confirms that more negative values of 
734 are found as ySb increases. The galaxy distribution in Fig. [3] 
is interpreted by making use of different dust-reddened stellar- 
population models. In particular, for a 100 Myr-old, solar metal- 
licity, continuous star-formation model of Maraston (2005 ), we 
show models with increasing Eg-v for the attenuation law of 
Calzetti et al. (120001) . different clumpy shell models of Witt 
& Gordon (|20001l. and uniform screen models using different 
extinction curves obtained for the SMC, LMC (Gordon et al. 
12003 ), and the MW (Cardelli et al. [T989] l. The distribution of the 
data points is well constrained by these models, and evidence 
of extinction curves with a form other than that of the SMC 
(with no UV bump) is clear The evidence of non-SMC extinc- 
tion curves is insensitive to the choice of the stellar-population 
model. For instance, considering a Maraston model with a SFR 
e-folding time of 3 Gyr at an epoch of 3 Gyr after the beginning 
of the SF activity would shift the plotted model curves only by 
H-0.20 in j6b and +0.08 in 734. 

According to the models shown in Fig. [3] the presence of 
a moderate UV bump (as in the LMC extinction curve) is sug- 
gested. However, the identification of the extinction curve be- 
comes complex, if models implementing age-dependent extinc- 
tion scenarios are considered. In fact, as a given simple stellar 
population becomes less attenuated as time goes by (i.e., as the 
average age of the stars increases), extinction curves with more 
pronounced UV bumps (e.g., as for the MW) are allowed by 
the distribution of the data points (see NPPS 'for a detailed dis- 
cussion). This scenario can be implemented by, e.g., realistic 
dust/star distributions similar to those of nearby disc galaxies 
(e.g., Silva et al. 1998, Pierini et al. '2004bi However, it cannot 
explain the existence of spectra with extremely negative values 
of 734. In this case, extraplanar dust must be invoked (NPPS). 
Hence, time-independent screen models (without scattering) or 
dusty shells appear to describe dust attenuation well in the high- 
redshift galaxies of our sample, and, thus, the data distribution 
in Fig. [3] 

Besides the overall agreement, some slight but interesting 
differences between samples exist, as suggested by Fig.|5] which 
reproduces average values of ySb for A734 = 1 bins. In fact, mov- 
ing from the FDF galaxies at 1.5 < z < 2.5 ((z) = 2.21) to 
the GMASS galaxies in the same redshift range ((z) = 1-89), 
and, then, from the FDF galaxies at 1 < z < 1.5 ((z) = 1.27) 
to the K20 galaxies in the same redshift range ((z) - 1.18), j6b 
tends to increase or the rest-frame UV continua become redder 
Comparison of values obtained for the different samples shows 
that the significance of this trend ranges between 1 and 3 cr. 
The most reliable values are found for the largest redshift differ- 
ences. In principle, differences in the selection criteria, average 
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Fig. 6. Comparison of composite spectra of 1.5 < z < 2.5 
GMASS galaxies with 734 > -2 (blue) and 734 < -2 (red), 
respectively. The lower panel gives the ratio of both composites, 
normalised at 2400 - 2570 A. Deviations from a Calzetti law 
in the UV-bump range are illustrated by the overplotted best-fit 
Calzetti attenuation curve with Eb-v - 0.08. 

luminosities, and redshifts (partly driven by the small number 
statistics available for individual samples) or systematic error^ 
could play a role. For the GMASS galaxies at 1.5 < z < 2.5, 
we checked that the variation in ySb is negligible (0.06 + 0.09) for 
objects brighter or fainter than R - 24. Therefore, the presence 
of redder UV continua as redshift decreases may have a physi- 
cal origin, i.e., ageing of the overall stellar population and/or an 
increase in the overall ISM opacity. A change in the dust com- 
position and/or distribution is possible as well. 

Interestingly, the sample of 24 local starburst galaxies ob- 
served by lUE and selected by INP05I exhibits relatively bluer 
UV continua than our sample galaxies: (/3b) - -1.79 + 0.11 
for (734) = 0.27 + 0.32. Suitable models of a MW-like galaxy 
that involve stellar population synthesis and radiative transfer 
for a disc geometry (see Pierini et al. I2003I I predict redder UV 
continua than starbursts in the local Universe. According to the 
models used by NP05 ageing of the stellar populations would 
however produce only a moderate shift in ySb equal to 0.1-0.2. 
Unfortunately, a robust age-dating of the stellar populations of 
all our sample galaxies is not possible, since the redshifted 
Balmer and 4000 A breaks are detectable in optical spectra only 
for z < 1.5 (see lNPPSl l. and model fitting provides quite uncer- 
tain age estimates. 

4.2. Constraining tlie attenuation curve 

In analogy with INP05 1 and INPPSFI. we divide the GMASS sam- 
ple of galaxies at 1.5 < z < 2.5 into subsamples, using the same 

^ The FDF galaxies dominate the sample of galaxies beyond z = 2.3. 
Their measured values of ySb can be affected by systematic errors due 
to the presence of strong OH residuals in the red optical wavelength 
region, where the redshifted /?(, was measured. 

The lNPPSl sample of galaxies at 1 < z < 1.5 is divided into a blue 
(fib < -1.5) and two red subsamples differing in 734. 
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Sample" 


Selection 








Ee- 


-V 


xo L"m '] 


^0 [A] 


7 \pm '] 


FWHM [A] 


(a) gm_wbb 


GMASS, 1.5<z<2.5, 7i4>- 


-2, /3b 


< 


-1.5 


0.164 ± 


0.069 


4.493 ± 


0.024 


2226 + 12 


0.498 ± 


0.113 


247 ± 


56 


(b) gm_wbr 


GMASS, 1.5 <z < 2.5, 7^ > - 


-2, /3b 


> 


-1.5 


0.286 ± 


0.071 


4.558 ± 


0.024 


2194 ± 11 


0.682 ± 


0.052 


328 ± 


25 


(c) gm_sb 


GMASS, 1.5 <z< 2.5,7^4 <- 


-2 






0.296 ± 


0.069 


4.572 ± 


0.010 


2187 ±05 


n -t- 

U. J 1 J ± 


U.UZl 


97A -t- 
Z /M- It 


1 n 

lU 


(d) n07_sb 


FDF/K20, 1.0 <z< 1.5,734 < 


-2 






0.344 ± 


0.078 


4.570 ± 


0.019 


2188 ±09 


0.583 ± 


0.083 


279 ± 


38 


(e) LMC2JOW7 


LMC 2 supershell (7 < 1) 








0.188 ± 


0.009 


4.540 ± 


0.053 


2203 ± 26 


0.684 ± 


0.042 


332 ± 


21 


Sample" 


Selection 








C3 [jim ^] 


^bump 


Abiimp [jjm-'] 










(a) gm_wbb 


GMASS, 1.5 <z< 2.5, 734 > - 


-2, /3b 


< 


-1.5 


0.116± 


0.055 


0.467 ± 


0.267 


0.366 ±0.192 










(b) gm_wbr 


GMASS, 1.5 <z< 2.5,734 >- 


-2, /3b 


> 


-1.5 


0.262 ± 


0.040 


0.562 ± 


0.106 


0.602 ±0.104 










(c) gm_sb 


GMASS, 1.5 <z< 2.5,734 <- 


-2 






0.305 ± 


0.020 


0.928 ± 


0.077 


0.835 ± 0.062 










(d) n07_sb 


FDF/K20, 1.0 <z< 1.5,734 < 


-2 






0.359 ± 


0.121 


1.056 ± 


0.186 


0.968 ±0.251 










(e) LMC2JOW7 


LMC 2 supershell (7 < 1) 








0.890 ± 


0.128 


1.903 ± 


0.320 


2.044 ± 0.320 











" For the samples at high redshift, the given errors are based on the uncertainties in the continuum of the mean spectrum without the 2175 A 
feature, Eb-v, and the FM-like parameter fit. For the LMC, they represent the sample variance. 



Table 1. Average properties of the three GMASS subsamples at 
1.5 <z < 2.5 



Par. 


734 > -2 


734 > 


-2 


734 < -2 




/3b < -1-5 


ySb > - 


-1.5 




A' 


29 


23 


16 


z 


1.86 ± 0.05 


1.85 ± 


0.05 


1.91 ±0.06 


R 


23.74 ± 0.09 


24.20 i 


:0.10 


23.99 ±0.13 




21.09 ±0.11 


20.74 4 


:0.15 


20.60 ±0.13 


734 


-0.20 ±0.16 


-0.62 H 


:0.19 


-3.02 ±0.19 


/Jb 


-1.83 ±0.04 


-1.24 H 


:0.04 


-1.18 ±0.10 


logLisoo [W/A] 


33.95 ± 0.04 


33.68 4 


:0.03 


33.82 ±0.07 


logL2500 [W/A] 


33.61 ±0.04 


33.46 4 


:0.04 


33.60 ± 0.06 


Eb-v 


0.133 ±0.020 


0.244 ± 


0.022 


0.304 ± 0.020 


logLboi [Lq] 


11.37 ±0.08 


11.564 


:0.07 


11.85 ±0.04 


logSFR[Me/yr] 


1.53 ± 0.09 


1.70 ± 


0.08 


2.05 ± 0.04 


logM, [Mo] 


9.94 ± 0.05 


10.12 4 


:0.08 


10.15 ±0.07 


log0 [Gyr"'] 


0.52 ± 0.08 


0.51 ± 


0.11 


0.85 ± 0.07 


re [kpc]" 


2.84 ± 0.48 


2.65 ± 


0.37 


2.52 ± 0.38 


"scr 


1.56 ±0.32 


1.79 ± 


0.44 


0.77 ±0.15 


b/a 


0.46 ± 0.03 


0.54 ± 


0.04 


0.46 ± 0.04 


CcAS 


2.57 ± 0.07 


2.63 ± 


0.07 


2.55 ± 0.08 


^CAS 


0.40 ± 0.03 


0.42 ± 


0.03 


0.43 ± 0.03 


/24,.m [/^Jy] 


37.3 ± 14.8 


39.1 ± 


10.4 


63.3 ±9.9 


logLg/jm [Lq] 


10.28 ± 0.07 


10.35 4 


:0.09 


10.73 ± 0.07 


logLiR [Lq] 


11.11 ±0.09 


11.21 4 


:0.11 


11.66 ±0.08 


log(LiR/Luv) 


0.50 ±0.10 


0.89 ± 


0.13 


1.25 ±0.12 



" from log Te 



thresholds of /?b - -1-5 and 734 = -2. In particular, we consider 
one GMASS subsample with 734 < -2, and two with 734 > -2 
that differ in terms of /3b. For 1.5 < z < 2.5, 16 GMASS galaxies 
are SUBGs, whereas 23 and 29 GMASS galaxies have 734 > -2 
plus red or blue UV continua, respectively. Average values of 734 
and ySb are listed in Table[T]for these GMASS subsamples. 

Figure|6]shows a comparison of the composite spectra of the 
GMASS galaxies at L5 < z < 2.5 exhibiting either weak evi- 
dence (at most) of a broad absorption feature centred on 2175 A 
in their spectra (734 > -2), or a secure one (734 < -2). The pres- 
ence of this feature in the second (more reddened) composite 
spectrum is evident by eye. 



4.2.1 . The shape of the UV bump 

Lorentzian-like fits are applied to the stacked rest-frame UV 
continua obtained from galaxies at L5 < z < 2.5 in the three 
GMASS subsamples, and from FDF and K20 galaxies at LO < 
z < L5 with 734 < -fl as described in Sect. 14.11 The result- 
ing FM-like parameters describing the morphology of the ob- 
served UV bump — i.e., centre xq (or Ao), 7, and amplitude cj — 
and their errors are listed in Table |2l with the maximum height 
£bump and area Abump of the observed UV bump, and the mean 
Eb-v of the best-fit Calzetti et al. (2000) attenuation law used as 
a baseline. We note that the values of Eg-v obtained from this 
Lorentzian-like fitting are consistent with those obtained from 
the SED fitting of GMASS galaxies (see Sect.lTTTi within the er- 
rors. However, we acknowledge that differences in Eb-v become 
higher for red galaxies with some signature of the UV bump in 
their spectra. 

Figure |7] shows the quality of the Lorentzian-like best-fit 
functions of the observed (residual) UV bumps obtained for the 
four subsamples. We note that radiative transfer effects only 
smooth out the UV bump in disc configurations (cf. Ferrara et 
al. 119991 Pierini et al. 2004b). Furthermore, no significant vari- 
ation in the width of the observed UV bump is found for the 
different radiative transfer models of Witt & Gordon (2000) for 
a given dust composition. Thus, we can assume that the residual 
UV bumps determined from the composite spectra of our galax- 
ies are closely related to those present in the extinction curve 
characterising their dust. As a first result, Fig.[8]indicates a com- 
fortably good correlation between the maximum height of the 
UV bump /ibump (dependent on the FWHM and amplitude) and 
the independently measured 734. Since the increase in Eg-v with 
decreasing 734 in our sample (see Table|2]i is obviously unable to 
remove this correlation, the increasingly negative values of 734 
are indeed associated with more pronounced UV bumps. 

The central wavelengths and widths of the residual UV 
bumps obtained for GMASS SUBGs at 1.5 < z < 2.5 and FDF 
and K20 SUBGs at 1.0 < z < 1.5 are almost identical. They 
have counterparts among those measured along different sight- 



^ Their composite spectrum is cut at 1800 A, since only a few spec- 
tra extend shortwards. Possible effects of the UV bump on the 1800- 
1900 A region have been considered. 
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Fig. 7. Best Lorentzian fits of the bump component of the av- 
erage extinction curves of the three GMASS subsamples at 
1 .5 < z < 2.5 (a - c) and a subsample of FDF and K20 galaxies at 
1.0 < z < 1.5 with < -2 (d). Moreover, the average 2175 A 
bump of a sample of measurements in the LMC 2 supershell re- 
gion close to 30 Dor showing low y is presented (e). More details 
about the samples and the fits can be found in Table |2] 

lines towards the MW and LMC0 (see Fig.|9]l. In particular, the 
most suitable local counterpart is the supershell region LMC 2 
close to 30 Dor, for which the average UV bump from different 
sightlines (excluding two outliers) is centred on 2203 + 26 A and 
has a FWHM equal to 332 + 21 A (cf. Gordon et al. |2003] ). For 
instance, the central wavelength and FWHM of the residual UV 
bump obtained for the composite spectrum of GMASS SUBGs 
are equal to 2187 + 5 A and 274 + 10 A, respectively. In general, 
the observed UV bumps of our high-redshift galaxies are rela- 
tively narrow, since they exhibit only about 60% of the typical 
widths measured for sightlines towards the MW and LMC. 

In addition, we note that the two composite spectra obtained 
for either the blue or red GMASS galaxies at 1.5 < z < 2.5 with 
no strong evidence of a UV bump in their individual spectra (i.e., 
with > -2) nevertheless exhibit a weak one. In either case, 
the central wavelength and width of the residual UV bump are 
not dissimilar to those obtained for galaxies with 734 < -2. This 



^ The average values of the FM parameters for the MW (Fitzpatrick 
& Massa ,2007 J and LMC (Gordon et al. [2003] > are: 2 180 A (central 
wavelength), 440 A (FWHM) or 0.93 ^m"' (7), and 3.0 or 2.7 (cj), re- 
spectively. 
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Fig. 8. Amplitude of the UV bump in the extinction curve ^bump 
(see Sect. I3.2l i versus 734 indicating the strength of the 2175 A 
feature in the spectrum for the three subsamples of GMASS 
galaxies at 1 .5 < z < 2.5 and the'NPPS'subsample of 1 < z < L5 
galaxies with 734. Labels are the same as in Fig. [T] and Table |2] 
but for an open circle, which marks the position of an extinction 
curve without any UV bump. The corresponding 734 is given for 
the stellar population parameters described in Fig. [3] 

tells us that some traceable dust similar to that in the supershell 
region LMC 2 seems to be present also in a non-negligible frac- 
tion of the blue and red GMASS non-SUBGs at 1 .5 < z < 2.5. 

4.2.2. The shape of the attenuation curve 

Our previous results imply that, at least for part of the high- 
redshift galaxy population, the functional form of the attenuation 
curve differs significantly from a Calzetti law. This is demon- 
strated by Fig.|6] which shows the best-fit Calzetti law for the ra- 
tio of the composite spectrum for 734 > -2 to that for 734 < -2. 
Neglecting possible differences in the age distribution of the stel- 
lar populations of both composite spectra, the resulting spectrum 
in the lower panel of Fig. |6] should reflect only a difference in 
reddening, if the Calzetti law holds universally. As expected, 
the best-fit Calzetti law with Eb-y - 0.08 fails to reproduce 
the UV bump region. However, the reliable wavelength range 
of the composites (from 1450 to 2900 A) is too narrow to al- 
low a significant deviation in the slope of the attenuation law to 
be detected. Moreover, this simple test does not provide any con- 
straints on the reliability of the wavelength-dependent amount of 
attenuation predicted by the Calzetti law. Investigating the shape 
of the attenuation curve across a wide wavelength range would 
require a thorough study of a much larger data set and probably 
more sophisticated modelling. This is beyond the scope of the 
present work and is defeiTed to a future paper. 

For now, we propose to modify the Calzetti et al. ( |2000| l law 
k{X) by adding a Lorentzian-like UV bump C3Dvi,,y(/l) (see Eq.[T] 
in Sect. 13. 2t and derive the wavelength-dependent attenuation in 
magnitudes by 

A{X) = Eb-v {k{A) + C3D,„,/^)). (2) 
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Fig. 9. FM-like parameters xq (the inverse central wavelenth 
of the UV bump) and y (linked to the inverse FWHM of the 
UV bump) for the three subsamples of GMASS galaxies at 
1.5 < z < 2.5 and the INPPSi subsample of 1 < z < 1.5 galax- 
ies with < -2 (filled symbols). Labels are the same as in 
Fig.|2]and Table|2] Moreover, crosses mark average FM param- 
eters for the Milky Way (Fitzpatrick & Massa 120071 and differ- 
ent LMC environments (Gordon et al. I2003I ). The latter are the 
average LMC ('LMC_ave'), the entire LMC 2 supershell region 
('LMC2_ss'), and the LMC 2 supershell region restricted to stars 
with y < 1 /im"' ('LMC2Jowy'). The errors given are measure- 
ment uncertainties of the FM parameters for the corresponding 
mean extinction curves except for 'LMC2 Jowy', for which sam- 
ple mean errors are provided. 



Since we have already used the Calzetti law as a baseline for the 
derivation of the FM-like parameters (see Table|2]i, this approach 
allows us to fit the composite spectra of our subsamples. The 
Calzetti law appears to work quite well as an average attenua- 
tion law for large samples of high-redshift galaxies. Luminosity 
and SFR estimates derived by assuming the Calzetti law show 
good correspondence to values obtained by independent meth- 
ods (e g., Reddy & Steidel ,2004, Erb et al. |2006bt Daddi et al. 
I2007bl ). Therefore, we assume that our modified Calzetti law 
describes, to zeroth order, on average dust attenuation in high- 
redshift star-forming galaxies with weak to moderate UV bumps. 
Characteristic values of xq, y, and C3 are given in Table|2l 

On the other hand, one may expect that the a priori use of 
a pure Calzetti law for fitting SEDs of any high-redshift galaxy 
introduces systematic errors whenever the extinction curve pro- 
duced by dust in a galaxy exhibits a significant UV bump. As 
a test, we constructed synthetic SEDs that are consistent with 
our data using Maraston ( 120051 ) models reddened by a modified 
Calzetti law, according to a case 'gm_sb' UV bump (see Table|2]i. 
The corresponding broad-band filter fluxes for different redshifts 
were then fitted by using the correct attenuation curve and a pure 
Calzetti law. The latter yielded about 5% higher masses and 
about 10% lower SFRs. For extreme but realistic cases, these 
discrepancies could be twice as large. The exact values depend 
on star-formation history, redshift (because of the differing fil- 
ter coverage of the UV bump), and, of course, dust obscuration. 




Fig. 10. Star-formation rate versus total stellar mass for the 
GMASS sample galaxies with 1.5 < z < 2.5. The three 734 
and ySb selected subsamples are marked by different symbols (see 
legend). Specific SFRs of 10, 1, and 0.1 Gyr ' are indicated by 
lines. The upper left corner of the diagram is not populated by 
galaxies because of the lower age limit of 100 Myr set for models 
with exponentially-declining SFRs for the SED fitting. 



The change in the rest-frame UV colours caused by the presence 
of a UV bump appears to be interpreted by a fitting routine that 
uses a standard Calzetti law as a signature of the older age of 
the stellar populations and/or a reduced dust reddening. For a 
typical high-redshift galaxy, the resulting deviations are almost 
negligible in view of the statistical and systematic uncertainties 
related to SED-fitting in general. In the following sections, we 
therefore discuss fitting results with the Hyperz code based on 
a standard Calzetti attenuation law to be consistent with other 
studies. Nevertheless, one should keep in mind that, for objects 
with the strongest 2175 A absorption features in their spectra, 
there could be noteworthy systematic deviations in the estimated 
galaxy properties. 

4.3. Extinction curve and global properties of galaxies 

4.3.1. Star-formation rates and stellar masses 

Table [T] lists average values of several global galaxy properties 
for the three GMASS subsamples defined in Sect. 14.11 In gen- 
eral, the monochromatic luminosities at 1500 and 2500 A (rest 
frame) of GMASS galaxies with a robustly identified 2175 A 
absorption feature in their spectra do not differ from the average 
ones of the entire GMASS sample. However, some differences 
exist in terms of the specific star-formation rate 0, which is the 
SFR per unit total stellar mass (derived from Hyperz, see 
Sect. Ol, as shown in Fig. [TO] In fact, GMASS SUBGs have 
<M*) = 1.4 X lO'^Mo and (SFR) ~ lOOMoyr-', so that their 
values of (p cluster around 7.1 Gyr"'. This is about twice the av- 
erage (p of GMASS galaxies for which the signature of a UV 
bump in their spectra is weak or absent ((0) = 3.3 Gyr"'), what- 
ever the observed UV slope. Interestingly, GMASS galaxies at 
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Fig. 11. Sersic index «.sei versus 734 for the three GMASS sub- 
samples at 1 .5 < z < 2.5 and three subsamples of FDF and K20 
galaxies at 1 < z < 1.5 discussed by .NPPS> The link between 
subsamples and symbols is illustrated in the legend. Open sym- 
bols indicate subsamples of galaxies with blue UV continua and 
weak 2175 A features. Errors represent sample mean errors. 



1.5 < z < 2.5 with > -2 and blue UV continua tend to be 
less massive systems with lower SFRs than corresponding non- 
SUBGs with red UV continua. 

A trend towards higher values of cp for SUBGs was also 
found bv lNPPSI We note that the GMASS galaxies at 1 < z < 
2.5 considered here are a factor of two less massive than the FDF 
galaxies at 2 < z < 2.5 in the lNPPSi sample, but almost as mas- 
sive as the FDF/K20 galaxies at 1 < z < 1 .5 in the NPPS sam- 
ple. Consistently, for the entire 4.5 jum-selected GMASS sample, 
which also comprises galaxies with photometry only and spec- 
troscopic redshifts from other projects in GOODS-S, the aver- 
age stellar mass drops from ~ 2.4 to ~ 1.0 x IO'^Mq moving 
from 2 < z < 2.5 to 1 < z < 1.5, if SFR > lOMo/yr (lower 
limit for our spectroscopic sample). Apart from different sam- 
ple selection criteria and different methods of estimating stellar 
masse the average sample masses probably decrease towards 
lower redshifts due to the increasing fraction of low-mass objects 
in flux-limited samples and a decrease in the number density of 
massive star-forming galaxies related to the "downsizing effect" 
(Cowie et al. 119961 Gavazzi & Scodeggio |1996l Gavazzi et al. 
[1996). 

4.3.2. Morphology 

Quantitative morphological studies of galaxies at z ~ 2 can be 
significantly affected by the relatively low surface brightness of, 
in particular, their outer regions, which can lead to erroneous 
morphological classifications. However, extremely deep ACS 



^ Comparison of stellar masses for galaxies in both samples suggests 
that the values used bv lNPPSl are systematically higher by about 0. 1 dex, 
which is probably in general due to the different assumptions about the 
SFH (see also Drory et al. l2005t . 



exposures are available for GOODS-South (see Sect. l2.1b . which 
allow us to analyse the rest-frame UV/U-band surface bright- 
ness distribution of our sample galaxies at least in an approxi- 
mate way. Table [T] shows the results for the effective radius r^, 
Sersic index Wser, elongation b/a, and the CAS parameters con- 
centration C and asymmetry A (see Cassata et al. 12005 1 for the 
definitions used) for our three GMASS samples at 1 .5 < z < 2.5. 

In general, the derived morphological parameters do not de- 
pend significantly on the UV continuum properties. A remark- 
able exception is the average Sersic index for GMASS SUBGs 
("ser) = 0.77 + 0.15, which is conspicuously lower than the val- 
ues 1.56 + 0.32 and 1.79 + 0.44 for coeval non-SUBGs with 
blue and red UV continua, respectively. Consequently, SUBGs 
appear to show radial distributions of the young stellar popu- 
lations that are on average shallower than an exponential pro- 
file. On the other hand, half of the SUBGs have the minimum 
value allowed for the fit Wser = 0.5, which challenges the re- 
liability of the fits. Moreover, there is no significant difference 
between our three GMASS subsamples in terms of the concen- 
tration C, which also depends on the radial surface brightness 
profile. However, this dependence is different, since C is related 
to the ratio of the apertures containing 80% and 20% of the total 
flux, while «ser describes the functional form of the radial light 
distribution. The concentration C is obviously less sensitive to 
the substructure of the radial profile than «sei, which appears to 
indicate irregular, very patchy rest-frame UV surface brightness 
profiles for low Sersic indices. A good argument for the relia- 
bility of the observed y34-related trend in Wser is the result that 
a similar relation is present in the intermediate-redshift NPPS] 
sample at 1 < z < 1.5 (see Fig. fTTT i. Even if the physical in- 
terpretation of the Hsei values is diflicult, the significance of the 
observed trend could in principle be used to search (in combina- 
tion with other tracers) for SUBG candidates without available 
deep UV spectroscopy. 

The use of the CAS parameters for galaxy classification (see 
Conselice et al. |2003l l shows that the GMASS galaxies at 1.5 < 
z < 2.5 appear to be a fair mixture of mergers/irregulars and 
disc galaxies, whatever 734. This distribution differs significantly 
from the one at 1 < z < 1.5, where most of the SUBGs appear 
to have late-type morphologies (see NPPS|. Since mergers are 
frequent in the NPPS sample of SUBGs at 2 < z < 2.5, this 
suggests that the distribution of the young stellar populations in 
galaxies with LMC-like dust becomes smoother as time goes by. 



4.3.3. Narrow spectroscopic features 

In the rest-frame UV domain mapped by the composite GMASS 
spectra, we can study 15 well-detected spectral lines from 
Si IV /1 1400 to MgU2852 (see Table |3j. The most prominent 
features are interstellar absorption lines and/or originate in stel- 
lar winds, and often show P Cygni profiles. Equivalent widths 
(EWs) and their uncertainties are reproduced in Fig. [12] as well 
as relative EWs with respect to the GMASS subsample of galax- 
ies at 1.5 < z < 2.5 with 734 > -2 and red UV continua. 
Comparison of the 12 absorption-dominated spectral lines of the 
GMASS blue and red subsamples of non-SUBGs reveals that the 
former type of galaxies has weaker features by only 1 1 (+5)% 
on average. In contrast, GMASS SUBGs exhibit stronger fea- 
tures, on average, by 34 (+5)%. On the other hand, the GMASS 
red subsample of galaxies with 734 > -2 tend to have stronger 
emission-dominated lines than the others by 20-30%. However, 
these emission lines have relatively uncertain EWs due to their 
weakness, even in composite spectra. 
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Table 3. Rest-frame equivalent widths in A of prominent lines in 
the composites of the three GMASS subsamples at 1.5 < z < 2.5 



Line Type" 


734 > 


_2 


/34 ^ ^ 


734 < ^ 






Pb ^ 


1.5 






Ol 1 V /t 14UU 


is 


+3.18 ± 


0.50 


+2.90 ± 0.63 


+4.06 ± 0.75 


o V /( ljUz 


s 


+0.55 ± 


0.09 


+0.60 ± 0.07 


+0.83+0.11 


c;tt 1 ICO/; 


i 


+ 1.97 ± 


0.12 


+ 1.95 + 0.17 


+2.73+0.28 


U 1 V /t 1 jjU 


si 


+4.48 ± 


0.36 


+4.52 + 0.54 


+5.65+0.32 


FelU1608 


i 


+ 1.36 ± 


0.15 


+ 1.72 + 0.12 


+ 1.72 + 0.27 


He lU 1640 


sn 


-1.14 ± 


0.27 


-1.23+0.35 


-1.29 + 0.28 


AlIU 1670 


i 


+ 1.43 ± 


0.08 


+ 1.53 ± 0.14 


+2.04 ± 0.18 


AlllU 1860 


i 


+2.61 ± 


0.20 


+3.11+0.22 


+3.66 + 0.20 


C 111] /1 1907 


n 


-1.05 ± 


0.14 


-1.48 + 0.14 


-1.16 + 0.22 


C 11] X ll^lb 


n 


-1.03 ± 


0.22 


-1.63 + 0.21 


-0.92 + 0.21 


FelU2344 


i 


+ 1.97 ± 


0.35 


+2.07 + 0.19 


+2.95 ± 0.26 


FelU2380 


i 


+3.16 ± 


0.32 


+3.67 ± 0.24 


+4.83 ± 0.29 


FelU2600 


i 


+4.79 ± 


0.48 


+5.15 + 0.62 


+5.85+0.68 


Mg 11 /1 2800 


is 


+4.01 ± 


0.75 


+4.02 ± 0.73 


+6.59 + 0.76 


MgU2852 


s 


+0.30 ± 


0.19 


+0.83 ± 0.25 


+ 1.30 + 0.34 


" Line origin 


s = 


stellar (photosphere and/or wind), 


i = interstellar 



absorption, n = nebular emission 



Hence, we confirm that galaxies with secure detections of 
a broad absorption feature centred on 2175 A in their spec- 
tra exhibit stronger low-ionisation interstellar absorption lines 
than those with a probable SMC-like extinction curve (INP05I I. 
This was interpreted as being due to a higher covering fraction 
of dusty neutral clouds in the former galaxies, which enables 
the carriers of the UV bump to survive. In this case, the low- 
ionisation interstellar lines trace the distribution of neutral gas 
clouds along the observer's sightline towards the hot stars pro- 
ducing the illuminating UV continuum radiation (Shapley et al. 
12003 ). An alternative explanation for the presence of stronger in- 
terstellar absorption lines is a broader velocity distribution (see 
Heckman etaL [T998] l. 

In addition. Fig. [12] reveals that, in galaxies with a signa- 
ture of a UV bump in their spectra, stronger absorption is in- 
ferred from all types of absorption lines, regardless of the wave- 
length and ionisation level. Hence, the higher absorption effi- 
ciency must be explained also for line components originating 
in stellar photospheres and winds. An enhanced metallicity is a 
likely interpretation. Since the strengths of wind features depend 
on metallicity (e.g., Leitherer et al. 1200 II IVIehlert et al. ,2002j , a 
higher metallicity could explain (at least in part) the significant 
absorption in the CIV /1 1550 A feature of SUBGs (cf^POS]!. 
At least for the red galaxies in our sample, this interpretation is 
supported by the coiTesponding presence of UV nebular emis- 
sion lines of lower strength (see Heckman et al. ,1998i Noll et al. 
I2004I I in the spectra of the same 734 < -2 galaxies. 

Unfortunately, a determination of metallicity for the com- 
posite spectrum of each GMASS subsample considered here is 
not possible. Halliday et al. (120081 ) were able to determine a 
stellar metallicity equal to ~ 0.3 Zq, based on the weak stellar 
Fe III /1 1978 A index (Rix et al. 120041 ). for a composite spec- 
trum combining 75 galaxy spectra from the entire GMASS sur- 
vey. This iron-abundance, stellar metallicity is lower than the 
oxygen-abundance, gas-phase metallicity determined for UV- 
selected star-forming galaxies at z ~ 2 (Erb et al. 2006a), and 
metallicities determined for hot stars (Mehlert et al. 2002 , 2006 ) 
and the ISM (Savaglio et al. 120041 ) in galaxies at similar red- 
shifts, which exhibit values ^ 0.5 Zq. Halliday et al. (120081 ) con- 
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Fig. 12. Rest-frame equivalent widths of 15 spectral fines be- 
tween 1350 and 2900 A for the composite spectra of three 
GMASS subsamples at 1.5 < z < 2.5. Circles and lozenges 
mark blue and red galaxies with 734 > -2. Squares indicate 
galaxies with 734 < -2. The 1 cr eiTor bars include the variance 
between the individual spectra, the errors of the continuum level 
definition, and the statistical uncertainties of the line strengths. 
Upper panel: Absolute equivalent widths in A. Emission lines 
have negative EWs. Lower panel: Equivalent widths relative to 
the values for the red subsample with weak 2175 A absorption 
features. Uncertainties on the reference values are considered 
by the error bars of the other samples. The relative EWs of the 
absorption-dominated lines of each subsample are connected by 
dotted lines in order to guide the eye. 



eluded that a light-element overabundance was being established 
in the selected GMASS star-forming galaxies as they were form- 
ing at z ~ 2. The higher EW of the Mg I A 2852 A line measured 
in GMASS galaxies that probably host the carriers of the UV 
bump could be explained by a higher fraction of intermediate- 
age stars (cf. NPPS). Unfortunately, uncertainties prevent us 
from establishing any diff'erence in either the measured Mg-to- 
Fe EW ratio or the properties of the stellar populations (from 
our models) among our GMASS subsamples for 1.5 < z < 2.5. 
Nevertheless, we note that the stellar metallicities of our star- 
forming, intermediate-mass galaxies at 1.0 < z < 2.5 are not far 
from those of LMC stars. 
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Fig. 13. Rest-frame luminosity at about 8 jim L%^^ versus 734 for 
the three GMASS subsamples at 1.5 < z < 2.5 (see legend for 
the symbols). Arrows indicate objects with 3 cr upper flux limits 
only. 



4.4. Dust grain populations emitting at mid-IR wavelengtlis 

The carriers of the UV bump probably consist of organic car- 
bon and amorphous silicates (Bradley et al. 2005 and refer- 
ences therein), and seem to populate different environments of 
their host galaxy (see, e. g., Wh ittet |2003 ; Gordon et al. 120031 
Fitzpatrick l2004l Clavton l2004l and references therein). Another 
component of the dusty ISM, which is ubiquitous and also con- 
sists of carbonaceous material, is represented by the large poly- 
cyclic aromatic hydrocarbon (PAH) molecules (see Peelers et al. 
|2004| and references therein). Far-UV pumped, vibrational emis- 
sion of PAHs is generally held responsible for the well-known 
emission features at 3.3, 6.2, 7.7, 8.6, and 1 1.2yum, which domi- 
nate the mid-IR spectra of many objects with associated dust and 
gas (Leger & Puget 1984; Al laman dola et al. 1985 1989; Puget 
& Leger 1989; Tielens et al. [19991 12000). Unfortunately, there 
is no mid-IR spectroscopy for the GMASS sample, so that the 
physical conditions in the global environment where PAHs and 
the caiTiers of the UV bump exist cannot be inferred either quan- 
titatively or qualitatively (from the comparison with the spectra 
of nearby objects). However, we can at least use the emission in 
the broad-band MIPS 24 yum filter to estimate the aUeged PAH 
emission around ~ 8 fim (rest frame). 

Observations in the MIPS 24 //m filter are available for the 68 
GMASS galaxies under study: fluxes range from 6 to 410yuJy, 14 
galaxies being 2 cr detections, and an additional four having only 
3 cr upper limits. Interestingly, the four galaxies without detec- 
tion exhibit 734 > -2 and /3h < - 1 .5. The MIPS 24 fim fluxes are 
converted into rest-frame luminosities at ~ 8 yum, where Lg/zm is 
calculated from the filter-averaged vfv (cf. Daddi et al. |2007a| ). 
Relatively small cosmological ^-corrections are obtained from 
the Dale & Helou (2002) IR templates, which are parameterised 
by the slope of the power-law distribution of dust mass over heat- 
ing intensity a0. Of the GMASS UV-luminous, intermediate- 
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Fig. 14. Ratio of the IR luminosity Lir and the UV luminos- 
ity Luv versus the UV-continuum reddening measure /3b for the 
three GMASS subsamples at 1.5 < z < 2.5 (see legend for the 
symbols). Lir is derived from the 24 fim fluxes by using the Dale 
& Helou (I2OO2I 1 models, while Luv is derived from SED-fitting 
in the wavelength range below 3000 A. Galaxies that show par- 
ticularly high luminosity ratios in comparison to other galax- 
ies with similar /3h are separated by a dashed line. The mid- 
IR-excess objects of Daddi et al. ( |2007al ) are indicated by big 
open squares. X-ray sources in the catalogue of Luo et al. (120081 ) 
are marked by big crosses. The X-ray source above the dotted 
line is not part of the Daddi et al. (2007a) sample. Theoretical 
curves are given for solar-metallicity Maraston (120051) models 
with either continuous SFR and age of lOOMyr (black) or with 
T - 3 Gyr e-folding time SFR and age of 3 Gyr (magenta or 
grey). Solid lines are for a pure Calzetti et al. ( 120001) law. Dashed 
lines show the results for a Calzetti law plus the 'gm_sb' UV 
bump in Table |2] 



mass galaxies at 1.5 < z < 2.5, those with a signature of the UV 
bump in their spectra appear to be more luminous (by ~ 0.4 dex) 
at ~ 8 fim (rest frame) than the other galaxies, even if only ob- 
jects with similar UV-continuum reddening (fib > -1 .5) and stel- 
lar mass (cf. Table [TJ are considered (see Fig. [T3] l. The 8 //m lu- 
minosities of the former galaxies are basically above 2x10'" Lq. 

For a clearer understanding of these results, we show the 
diagnostic diagram defined by the luminosity ratio Lir/Lu\Q 
and the UV-continuum slope /3b for the 68 GMASS galaxies 
at 1.5 < z < 2.5 under study in Fig. [14] Most of these galax- 
ies are either located along the locus of nearby starbursts (cf. 
Meurer et al. '1999) or populate the region below this locus, 
where nearby "normal" star-forming galaxies lie (e.g., Kong et 
al. [20041 Buat et al. [20051 Panuzzo et al. [2007] l. According to the 



Templates are chosen according to the mean relation between a (or 
f^lfl ) and the IR luminosity derived by Marcillac et al. ( i2006» for 



local IR-luminous galaxies. We note that the Dale & Helou ( I2002t IR 
templates exhibit only small differences at 8;um (rest frame), so that 
small uncertainties are propagated to Lg^m- 

' LiR is derived from the measured MIPS 24yum fluxes by the adop- 
tion of the Dale & Helou ( I2002t templates and the calibration of 
Marcillac et al. J2006b . Luv is measured in the best-fit SED of Hyperz 
at wavelengths below 3000 A. 
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interpretations given by the previous authors, the scatter of the 
GMASS SUBGs across the locus of nearby starbursts (the so- 
called "Meurer law", see Meurer et al. 1 19991 1 is consistent with 
the simple dust-screen model described in Sect. 14. 11 At the same 
time, the location of (red) galaxies with 734 > -2 in Fig. [14] 
could indicate that the observed UV bumps are weakened due 
to more complex dust/star distributions (see lNPPSj for a detailed 
discussion). 

The a priori exclusion of optically identified AGN, still al- 
lows the presence of continuum emission at mid-IR wavelengths 
caused by deeply obscured non-thermal sources. However, only 
two galaxies among the GMASS galaxies plotted in Fig.[l4lhave 
been detected in X-rays: none of them exhibits a significant UV 
bump. In fact, a non-detection in X-rays does not exclude a hid- 
den AGN being present in our sample, as the X-ray stacking 
analysis of Daddi et al. ( I2007bl l indicates. Indeed, hidden AGN 
candidates identified as "IR-excess objects" (see Daddi et al. 
I2007aj are present among GMASS galaxies at 1.5 < z < 2.5 
with < -2 (25%) or with 734 > -2 (10%13- However, the 
moderate Lir/Luv for most of the GMASS galaxies in Fig. fT4l 
suggest that the measured MIPS 24 fim fluxes are mainly associ- 
ated with PAH emission and star formation. Moreover, our three 
GMASS subsamples do not show suspicious values and signifi- 
cant differences in IRAC colours that can discriminate between 
galaxies dominated by star formation and by AGN in the rest- 
frame near-IR (see Lacy et al. 2004; Stern et al. 2005; Pope et al. 
I2OO8I 1. For example, SUBGs and non-SUBGs show on average 
^8.0/5 4.5 - 0.8 and no SUBG was found to be located above the 
Po pe et a l. ©OS) AGN criterion of 2.0. Finally, the Teplitz et 
al. (I2007I I Spitzer IRS spectrum (source 2 in their Fig. 5) of the 
galaxy in the INPPSi sample with the strongest 2175 A feature, 
a K20 galaxy at z = 1 .097, indicates a typical, star-formation- 
dominated, luminous infrared galaxy. Unfortunately, this con- 
vincing check cannot be performed for other sample galaxies 
because of the lack of published mid-IR spectroscopy. 

Additional support against an AGN-related origin of the ob- 
served mid-IR emission comes from the stacking analysis of 
MIPS fluxes in the 70 /im filter for the three GMASS subsam- 
ples (D. Elbaz & B. Magnelli, priv.comm.). As a result, only 
3 cr upper limits between 600 and 800 //Jy are obtained. This 
excludes the suitability of Dale & Helou ( 2002 J IR SEDs with 
power-law slopes of the dust mass distribution over heating in- 
tensity a < 1 .2, which peak at wavelengths < 63 /^m in fy and are 
more extreme than SEDs of typical star-forming galaxies (Dale 
& Helou |2002] l. The typical values of a that we determine from 
the MIPS 24yum fluxes and the calibration of Marcillac et al. 
( |2006] | are about 1.7 - 1.8. 

We can therefore safely conclude that the bulk of the ob- 
served MIPS 24yum flux of the GMASS galaxies under study 
comes from PAH emission, this emission being higher for ob- 
jects with relatively strong 2175 A absorption features in their 
spectra. This enhanced PAH emission is probably due to at least 
one of the following reasons: 1) a higher abundance of PAHs or 
a larger amount of dust; 2) a more eflicient attenuation by dust as 
a result of the dust/star distribution (see lNPPSj for a discussion); 
3) a stronger illuminating radiation field, as a consequence of a 
larger UV throughput or a reduced dilution (i.e., increased close- 
ness of dust to illuminating stars). All of these independent inter- 
pretations are supported by our data to some level. In particular, 
the relation between the 8 fim emission and the SFR or (p (see 
TablellJ suggests a link between star-formation intensity and the 



One of the latter objects is actually part of the GOODS-S X-ray 
catalogue of Luo et al. 120088 (see Fig.ll4b. 



amount of dust in the galaxy (cf. Wang & Heckman ll996t Vijh 
etal. |2003l l. 



5. Discussion 

We have investigated spectral evidence of the broad absorp- 
tion feature centred on 2175 A, following the method devel- 
oped by NPO5] The INPP S main sample of 88 actively star- 
forming, intermediate-mass (i.e., with average stellar masses 
of 10"' - 10" M0) galaxies at 1 < z < 2.5 with appropri- 
ate spectra was enlarged by considering 68 GMASS galaxies 
at 1.5 < z < 2.5. The gain consists not only in the larger size 
and the more homogeneous redshift coverage of the new sample. 
The additional GMASS galaxies have particularly high S/N op- 
tical spectra that bracket the redshifted wavelength region across 
the UV bump, so that a Lorentzian-like fitting of this broad fea- 
ture can be performed. This allows us to study the shape of the 
UV bump and the link between its spectral signature (i.e., the pa- 
rameter 734 introduced by'NP05) and the properties (e.g., central 
wavelength and FWHM) of the fitted UV bumps. Furthermore, 
the GMASS sample is large enough to enable subsampling: in 
composite spectra of sufficiently high S/N, EWs can be mea- 
sured more robustly for many narrower spectral features than 
in previous studies. These absorption/emission lines originate in 
very different local environments, such as stellar photospheres, 
wind regions, nebular regions, and the ISM. Finally, available 
MIPS 24 yum observations allow us to estimate the total amount 
of emission by dust at IR wavelengths and constrain the opac- 
ity characterising different GMASS subsamples. Hereafter, we 
discuss our main results in terms of galaxy structure, and the 
evolution of stellar populations and the ISM. 

5. 1 . Differing extinction curves in high-redshift galaxies 

We have produced robust evidence that at least 30% of the UV- 
luminous, intermediate-mass galaxies at 1 < z < 2.5 contain 
dust that produces an extinction curve with a broad absorption 
excess centred on 2175 A that differs from the standard SMC 
curve at UV wavelength^ We refer to these objects, identified 
by 734 < -2, as significant L'V bump galaxies (SUBGs). With 
the exception of our previous investigations (NP05; NPPS), the 
presence of the caiTiers of the UV bump in statistical samples of 
galaxies has been limited so far to intervening Ca II and possibly 
Mg II absorbers at intermediate/high redshifts (Malhotra 119971 
Wang et al. 2004; Wild & Hewett 2005). Our results imply that 
a non-negligible fraction of the observed actively star-forming 
galaxies at 1 < z < 2.5 akeady has a UV extinction curve that 
differs from that characterising LBGs at similar and higher red- 
shifts (cf. Vijh et al. l2003l) . UV extinction curves in-between the 
average SMC one (with no UV bump but a steep far-UV slope) 
and the average LMC one (with a moderate UV bump and a 
less steep far-UV slope) appear to be common, whereas the av- 
erage MW extinction curve (with a prominent UV bump and a 
less steep far-UV slope) does not. In particular, the shape of the 
UV bump infeiTed from the composite spectrum of SUBGs in 
the GMASS sample resembles that measured in the "supergiant" 
shell of ionised filaments LMC 2, close to 30 DoiPl 



' ' Sample selection criteria and survey sensitivity can hinder this fact, 
as suggested by INPPSI 

This supergiant shell probably results from the combined action 
of stellar winds and supernova explosions over a timescale of 10^ yr 
(Caulet et al. J982i) and a unique magnetic field geometry, consisting of 
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Both the fraction of actively star-forming galaxies at 1 < 
z < 2.5 with alleged UV bumps in their spectra, and the 
type of extinction curve (e.g., LMC versus MW average ones), 
depend on the dust/star distribution, which is generally un- 
constrained by existing observations of high-redshift galax- 
ies (NPPS). Extinction curves with a moderate UV bump are 
favoured if dust is mainly distributed above stars, as suggested 
for SUBGs, which are indicative of high covering fractions of 
young massive stars due to large EWs of strong interstellar ab- 
sorption lines (Fig. [T2] ). Moreover, the distribution of GMASS 
SUBGs in the ^Sb-i-iR/iuv plane (Fig. [T4l i is aligned mainly 
along the opacity-dependent Meurer curve (Meurer et al. 1 19991 1 
for starburst models with screen-like dust attenuation. Dust 
screens are expected if winds on a galactic scale or more lo- 
calised superwinds are present (see, e.g., Murray et al. 120051) . 
Galactic outflows are ubiquitous, particularly at high redshift 
(e.g., Weiner et al. |2009t see Heckman 2005 and Veilleux et al. 
120051 for reviews). The presence of galactic winds or localised 
superwinds in GMASS galaxies, for instance, is supported by 
the high values of the specific SFR (Fig.fTOb. Part of the expelled 
dust grains may settle in a potential minimum within the galaxy 
halo, as a result of simple dynamical arguments about the mo- 
tion of interstellar grains out of the galactic disc (Greenberg et al. 
[19871 Barsella et al. [19891 Ferrara et al. [19901 [TWTl Davies et al. 
11998 1). The presence of poloidal magnetic fields may help to con- 
fine expelled dust grains (Beck et al. 1994), and, thus, preserve 
the SHELL-like configuration for the dust/star distribution that 
can describe dust attenuation in nearby and high-redshift star- 
bursts (see Witt & Gordon [20001 Vijh et al. l2TO3l) . Increasing 
observational evidence points to the existence of organised mag- 
netic fields in and around galaxies without AGN activity at inter- 
mediate/high redshifts, which are at least of comparable strength 
to those observed in nearby normal galaxies (Bernet et al. 120081 
Wolfe et al. |2008l l. 

A higher fraction of objects can host the carriers of the UV 
bump, if dust is distributed mainly in the plane of the galaxy (as, 
e.g., in the radiative transfer models presented in Silva et al. 119981 
and Pierini et al. |2004bj ). As shown bv lNPPSl a disc system with 
a MW extinction curve then exhibits -3 < 734 < 1 only. The 
fraction of dust above the mid-plane of a galaxy can decrease 
as the SFR per unit area of the galaxy decreases (cf. Heckman 
120051 ). Disc-dominated dust distributions become more likely as 
the cosmic SFR density decreases, i.e., towards redshifts pro- 
gressively lower than unity (Madau et al. 1 19961 ). On the other 
hand, extinction curves more similar to those of the LMC or MW 
are also possible, if the dust/star distribution is patchy on very 
large scales (see INPPSI ) or depends on the age of the individ- 
ual stellar populations (see Panuzzo et al. 2007). These dust/star 
distributions appear to explain the position of "normal" star- 
forming galaxies distinctly below the Meurer curve (Meurer et 
al. |1999t Kong et al. 2004, Buat et al. 2005 ), where many non- 
SUBGs in our sample are located (see Fig.[T4|. 

Finally, we note that the extinction properties inferred in this 
study and in our previous ones refer to a galaxy-averaged envi- 
ronment. It is well known from studies of the MW, LMC, and 
SMC that bulk dust properties vary on much smaller scales than 
a galaxy scale length (down to sub-pc). Hence, the carriers of the 
UV bump may exist in particular regions of galaxies that do not 
exhibit spectral evidence of the 2175 A feature. This may explain 
why a typical LMC extinction curve has so far been inferred only 
for the sightline towards GRB 070802 at z = 2.45 (Eli'asdottir et 



coherent fields on spatial scales of 42x24 pc to 104x83 pc with inferred 
projected field strengths of ~ 14 - 30yuG (Wisniewski et al. 12007) . 



al. 120081 ), whereas GRB-host galaxies are commonly assumed 
to host SMC-like dust. It is clear that the Stardust yields and 
the processes of creation and transformation of dust species and 
sizes in the general ISM play a fundamental role in explaining 
the variety of extinction curves in different galaxy environments. 
This is discussed in the next section. 

5.2. The carriers of the UV bump and other small-size dust 
components at high redshifts 

Some forty years after the discovery of the UV bump, pro- 
posals about the chemical composition of its carriers remain 
largely conjectural. If the observed properties of the UV bump 
strongly suggest that its carriers are in the small particle limit 
(i.e., a « A), the abundance argument limits the set of possible 
constituent elements to C, Mg, Si, and Fe. Graphitic carbon is 
by far the most widely discussed and accepted material among 
these candidates, since it is the only one to exhibit a broad excess 
in the absorption cross-section at UV wavelengths (see Papoular 
& Papoular 2009 for a recent discussion). However, as for any 
material, its properties must explain the variations in width of 
the UV bump, whilst conserving its peak wavelength (Whittet 
12003 0- The random, hydrogen-free assembly of microscopic 
sp^ carbon chips forming a macroscopically homogeneous and 
isotropic solid, which was proposed as a model carrier of the UV 
bump (Papoular & Papoular 2009) , is able to explain these prop- 
erties taking into account astrophysical conditions in the ISM. 
The most promising, alternative carbon-based materials are syn- 
thetic nanoparticles condensed in a hydrogen-rich atmosphere 
and polycyclic aromatic hydrocarbons (PAHs), which exhibit a 
rise in the absorption cross-section across the region in which 
the UV bump originates. In addition, oxygen-rich materials have 
been proposed as carriers of the UV bump. In particular, silicate- 
based materials appear the most probable candidates, although 
this may require unreasonable fine-tuning of the size distribution 
(Whittet 2003 ). If graphitic grains (at least partly) are responsi- 
ble for the UV bump, the main effect of the presence of PAHs 
(and probably silicates) will be then to increase the absorption 
cross section in particular at the wings of the UV bump, caus- 
ing the feature to become less prominent. A variable mixture of 
the three components could explain variations in the height and 
width of the UV bump among observed extinction curves. Thus, 
it is safe to assume that the caiTiers of the UV bump are a mix- 
ture of amorphous silicates and carbonaceous material, as indeed 
observed in solar-system interplanetary particles (Bradley et al. 
20050 

According to the successful dust model of Desert et al. 
( 119901 ), the far-UV non-linear rise in the extinction curve can 
be ascribed mostly to large molecules (radius a ~ 0.4 - 1 nm). 



The width of the UV bump possibly varies owing to different mech- 
anisms, including particle clustering, ice-mantle growth, compositional 
inhomogeneities, porosity variations, and surface effects. Our data do 
not give us a handle on any of these processes; however, we expect less 
particle clustering and mantle growth in carriers that are above the mid- 
plane of the galaxy (as in SUBGs), since the two phenomena take place 
in regions like molecular clouds, protostellar discs, etc. 

Carbon-rich and silicate materials are made in different circumstel- 
lar environments (see text later on). How and where the two compo- 
nents, originating in different production sites, were combined within 
the same dust particles remains unclear. Hence it is uncertain whether 
the history of solar-system particles can be regarded as typical of those 
in remote galaxies. Moreover, Bradley et al. 120051 could not distin- 
guish whether the UV bump is due to organic carbon or amorphous 
silicates (or both), nor could exclude other possible UV-bump carriers. 
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i.e., PAHs, whereas carbonaceous nanoparticles (a ~ I - 10 nm) 
called very small grains (VSGs) are responsible for the sharp ab- 
sorption excess at 2 1 75 A, and big grains (BGs, a ~ 1 - 1 00 nm) 
cause most of the extinction at UV wavelengths, but in a fea- 
tureless way and with decreasing efficiency from the near- to the 
far-UV domain. Furthermore, PAHs are the caiTiers of the aro- 
matic infrared bands between 3 and 17 fj.m, whereas VSGs are 
responsible for the continuum emission at mid-lR wavelengths. 
Differences in the far-UV slope and/or the strength of the UV 
bump are related to zeroth order differences in the fraction of big 
grains and in the VSG-to-PAH abundance ratio of the dust. The 
latter arise from both dust production, either in stellar outflows 
(i.e., Stardust) or in the general ISM, and dust processing in the 
general ISM (Whittet [2003l and references therein). Regardless 
of whether the carriers of the UV bump originate in stars, the 
main unknown is the abundance required to alter in the appro- 
priate way the shape of the average extinction curve of a galaxy. 

Silicates (mostly amorphous) are produced in outflows of 
oxygen-rich cold stars, such as O-rich asymptotic giant branch 
(AGB) stars and red giant and supergiant stars (Whittet |2003l 
and references therein). In contrast, amorphous carbonaceous 
grains and PAHs are produced in outflows of carbon-rich cold 
stars, i.e., AGB carbon star£3- In any event, intermediate-mass 
stars (1 < M < SM©) lose most of their mass during their 
AGB phase (e.g., Iben & Renzini 1983). Hence, they are ma- 
jor producers of dust grains, possibly including those responsi- 
ble for the UV bump. Besides the winds of AGB stars, the ex- 
panding envelopes of Type II supernovae (SNe) represent an O- 
rich environment where silicate dust can form (Whittet 2003). 
According to Zhukovska et al. (2008 ), SNe could also be im- 
portant for the production of carbon dust, mostly in crystalline 
form (i.e., graphitic-like). In any case, the importance of SNe for 
the dust production in galaxies is still highly uncertain (Bianchi 
& Schneider 2007 and references therein). AGB stars and their 
outflows appear with a short delay after the beginning of star 
formation, i.e., some ~ 30Myr, which is the lifetime of a 8 M© 
star. AGB stars more massive than ~ 3 Mq never become carbon 
stars (e.g., Renzini 2008 and references therein), can hence po- 
tentially produce silicate Stardust but not C-rich Stardust. For the 
latter grains to appear, one needs to wait for stars less massive 
than ~ 3 M© to reach the AGB, or ~ 300 Myr after the beginning 
of star formation. Moreover, low metallicity favours the appear- 
ance of carbon stars, because less carbon needs to be dredged up 
in AGB stars to reach C/O ratios of above unity, and turn them 
into carbon stars (Renzini & Voli |1981| ). In summary, the theory 
of stellar evolution foresees that silicates are more promptly re- 
leased than carbon-rich grains, and carbon grains are more abun- 
dantly produced in a low-metallicity environment. 

With the persistent uncertainty in the chemical nature of the 
UV-bump carriers, their stellar producers, and their possible re- 
processing in the general ISM, much of the following discussion 
is speculative. It is safe to say that the EWs of several emis- 
sion/absorption features of carbon and silicon imply that the 
chemical enrichment of these dust constituents has proceeded 
significantly in the SUBGs in our GMASS sample (see TableO. 
Additional evidence is provided by the emission at 8 yum (rest 
frame), if this emission is dominated by PAHs. In particular, 
emission features around 8 fim (rest frame) are associated with 



The amorphous carbon grains ejected by C-rich stars need to be at 
least partly graphitised by interstellar processes to produce the distinc- 
tive absorption feature at 2175 A (Whittet 2003 and references therein). 
This would take about 100 Myr if such a process is driven by UV radi- 
ation (see Sorrell |1990t Mennella et al. ri996ilT998] >. 



C-rich stars (Whittet [20031 and references therein). Thus, AGB 
outflows probably play a dominant role in the build-up of the car- 
riers of the UV bump in massive galaxies at high redshift (NPOSl 
NPPS ), as well as in the build-up of PAHs in local (Galliano 
2008 ) and high-redshift galaxies alike. Consistently, dust exclu- 
sively formed in the ejecta of core-collapse SNe and exposed to 
sputtering by the hot gas appears to produce extinction curves 
that do not exhibit a UV bump (see Bianchi et al. 2007; Nozawa 
et al. l2008l and references therein). 

The absence of the UV bump in the SMC challenges the idea 
that carbon grains represent its carriers, since the SMC is very 
rich in carbon stars and yet no UV bump is present there. The 
lack of a UV bump could also be caused by efficient grain de- 
struction (Gordon et al. 120031 Sloan et al. 120081 ) favoured by 
lower dust column densities and harder radiation fields due to 
the low metallicity. If carbonaceous carriers of the UV bump 
are more sensitive to harsh environments than silicate ones, then 
the UV bump strength may be lower in carbon-dust dominated 
low-metallicity environments than in high-metallicity galaxies, 
where the dominant population of silicates may also act as car- 
riers of the UV bump. The latter is possible for our sample of 
SUBGs (cf. Sect. |4.3.3t . On the other hand, these galaxies ap- 
pear to contain a large amount of PAHs (see Sect. 14.41 ). which 
are carbonaceous in nature. 

The destruction of refractory dust proceeds predominantly 
in the intermediate warm phase of the ISM (McKee 1989; Jones 
et al. 119941 ). This phase of the ISM has probably a large fill- 
ing factor in actively star-forming galaxies at high redshift. In 
the low-density phases of the ISM, where SN explosions gen- 
erally occur (at least in the local Universe), the relative abun- 
dance of PAHs with respect to VSGs is higher than in the high- 
density phases of the ISM. This was also found to be true for 
the diffuse, illuminated part of dense photodissociation regions 
(Compiegne et al. 2008 ), whose chemistry is controlled by the 
external UV radiation field. Thus, we can assume that, in diffuse 
(dense) regions of the ISM of UV-luminous, massive galaxies 
at high redshift, narrow (broad) UV bumps are observed as a 
consequence of the higher (lower) fraction of PAHs. This trend 
is observed in the MW (Fitzpatrick & Massa 119861 Cardelli 
& Clayton 119911 Valencic et al. 120041 ). although it is usually 
explained in terms of ice mantle growth in the dense regions 
(Mathis 119941 ). Unfortunately, a similarly clear trend has not 
been observed in the LMC (Misselt et al. 1999). On the other 
hand, different dust components may have different spatial dis- 
tributions on scales ranging from pc to kpc, and, thus, be exposed 
to different stresses caused by the UV radiation field and shocks. 
Unfortunately, the few existing computations of dust expulsion 
and erosion reach conflicting conclusions about the distances at 
which graphite and silicate grains are located and the equilib- 
rium distances reached by grains of different sizes (cf. Bianchi & 
FeiTara 2005, AguiiTe et al. I2001II200 1). Our analysis suggests 
that caiTiers of the UV bump populate regions at some distance 
from the mid-plane of a galaxy at least in SUBGs. They might 
have survived expulsion or been produced by erosion during ex- 
pulsion of larger carbonaceous grains. 

Whatever the carriers of the UV bump are, it is reasonable to 
expect that a higher fraction of them are hosted by more evolved 
galaxies at high redshift. If they are mostly silicates, galactic nu- 
cleosynthesis would be required to accumulate sufficient silicon 
and oxygen and enable them to condense onto grains in mas- 
sive AGB star winds. If they are mostly carbonaceous, stars of 
lower mass than ~ 3Mq have to reach the AGB (~ 300 Myr) 
and start producing C grains, their production rate peaking when 
~ 2 Mq stars reach the AGB, about 1 Gyr later. If the carriers 
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of the UV bump are not Stardust, they are produced in the gen- 
eral ISM on timescales of the order of 0.5 Gyr (Whittet 120031 
and references therein). The dust production timescales can be- 
come significantly shorter, if SNe of Type II (which could occur 
at very high redshift; see Liang & Li '2008; Nozawa et al. l2008l l 
or the slow winds of their precursors (red supergiants) play a 
role. In any case, the carriers of the UV bump have to spread 
and survive across and above the disc (possibly as high as 10 
kpc), where part of them (or their parent grains) are transported 
to speeds of 10 to 500km s"' by galactic winds and superwinds. 
Finally, if the presence of the carriers of the UV bump is ruled 
not only by their accumulation (in stars or in the ISM) but also 
by their survival (e.g., favoured by self-shielding; see NPOS,), 
some threshold in either dust-to-gas ratio or metal abundance 
may have been reached (see Gordon et al. 120031 1. 

6. Conclusions 

We have selected a total of 78 UV-luminous galaxies at 1 < z < 
2.5 with average stellar mass of about 10'° Mq from GMASS, a 
deep spectroscopic survey of IRAC 4.5 /zm-flux limited galaxies 
at z ~ 2 in the GOODS-S (Kurk et al. |2008aD . This sample ex- 
tends to 186 the number of actively star-forming, intermediate- 
mass (i.e., with stellar masses of 10'" - 10" Mq) galaxies at 
1 < z < 2.5 with high S/N optical spectra and multiwave- 
length photometry, selected from the FDF, K20, and GDDS sam- 
ples and investigated bv lNPPSI Its heterogeneous nature assures 
us that the total sample probes the population of UV-luminous, 
massive galaxies at these redshifts comprehensively, beyond the 
limitations of the selection criteria of individual samples. 

The properties of the UV extinction curve, i.e., far-UV slope 
and presence (or absence) of the broad absorption feature cen- 
tred on 2175 A ("UV bump"), were constrained by a parametric 
description of the spectra by mapping the rest-frame UV spec- 
tral energy distributions of the individual galaxies. We have in- 
terpreted the data by applying a suite of different models com- 
bining radiative transfer and stellar population synthesis. For the 
68 GMASS spectra at 1.5 < z < 2.5 in particular, additional 
Lorentzian-like "Drude" fitting of the UV bump was performed. 
The overall analysis placed constraints on the dust population of 
large molecules and small grains responsible for the extinction 
at UV wavelengths. 

Hereafter, we summarise the main results of our analysis and 
place them in the context of our previous results. 

- About 30% of the spectra of actively star-forming, 
intermediate-mass galaxies at 1 < z < 2.5 undoubtedly ex- 
hibit broad absorption features centred on ~ 2175 A, regard- 
less of redshift. We term such galaxies SUBGs (.significant 
UV bump galaxies). Taking into account their measured UV 
slopes, we conclude that the extinction curves characteris- 
ing these galaxies differ from those observed in the SMC, as 
well as from those characterising nearby starburst galaxies 
and typical LBGs at 2 < z < 4. 

- If attenuation at UV wavelengths is dominated by dust dis- 
tributed above the stellar body of a galaxy, as suggested at 
least for SUBGs, effective UV extinction curves in-between 
those typically observed in the SMC and LMC characterise 
UV-luminous, intermediate-mass galaxies at 1 < z < 2.5. 
This dust/star configuration originate from the action of 
galactic winds or localised superwinds. Stronger intrinsic 
UV bumps are possible for an enhanced mixing of dust 
and stars (probably in a disc), patchier dusty screens, age- 
dependent extinction effects, and significantly diverse re- 



gions in terms of star formation and dustiness within a 
galaxy. Since this cannot be excluded, especially for non- 
SUBGs (for which the ratio Lir/Luv is significantly lower 
than expected for dust-screen dominated galaxies), the frac- 
tion of high-redshift galaxies with LMC-type extinction 
curves can be higher than inferred from the apparent bump 
strength. 

- The UV bumps of the analysed UV-luminous high-redshift 
galaxies exhibit about 60% of the width typically found in 
the LMC and Milky Way. In the local Universe, similar 
widths are found only for some sightlines towards the su- 
pergiant shell of ionised filaments LMC 2 near 30 Dor. This 
suggests similar environments in terms of radiation field, 
dust-to-gas mass ratio, and population of dust grains and 
molecules. 

- The presence of multiple carriers of the UV bump (probably 
organic carbon and amorphous silicates) in a galaxy argues 
that AGB stars have already started to return large amounts 
of dust to the ISM. The required chemical enrichment is par- 
ticularly evident for GMASS SUBGs at 1.5 < z < 2.5, for 
which relatively large equivalent widths can be measured for 
several UV absorption features, mainly of interstellar origin, 
produced by Si, S, C, Al, Fe, and Mg. 

- SUBGs are more luminous at 8 fj.m (rest frame) than most of 
the other GMASS galaxies at 1.5 < z < 2.5, which, on the 
other hand, have lower average star-formation rates. This is 
circumstantial evidence that carbonaceous dust has already 
been synthesised in relatively large amounts at these red- 
shifts, since polycyclic aromatic hydrocarbons (which origi- 
nate in AGB stars) probably dominate the emission at 8 yum. 

- Large amounts of UV radiation and shocks must be present 
in our sample galaxies. This is particularly true for SUBGs, 
since their specific star-formation rates are higher. Survival 
of the small-size dust component in such harsh environments 
is helped by the higher injection rates and more efficient dust 
self-shielding. Consistently, UV-luminous, massive galaxies 
at 1 < z < 2.5 with LMC2-like dust suffer from higher at- 
tenuation at UV wavelengths and, according to our models, 
host larger amounts of dust. 

The frequent occurrence of significant UV bumps in spec- 
tra of high-redshift galaxies appears to imply that these galaxy 
SEDs could be reproduced more accurately by models with at- 
tenuation curves that differ from the Calzetti law, which does 
not account for the presence of a UV bump. The investigation 
of the attenuation curves of star-forming, high-redshift galaxies 
across a wider wavelength range is deferred to a future paper. 
However, the present study indicates that a Calzetti law plus a 
Lorentzian-like bump reproduces the spectra of SUBGs at least 
for the wavelength range 1450-2900 A. Accounting for the pres- 
ence of a UV bump in this way involves systematic but small 
changes in crucial properties such as the SFR in SUBGs. 
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